Drought trends projection under future climate change scenarios for Iran region

The study highlights the potential characteristics of droughts under future climate change scenarios. For this purpose, the changes in Standardized Precipitation Evapotranspiration Index (SPEI) under the A1B, A2, and B1 climate change scenarios in Iran were assessed. The daily weather data of 30 synoptic stations from 1992 to 2010 were analyzed. The HadCM3 statistical model in the LARS-WG was used to predict the future weather conditions between 2011 and 2112, for three 34-year periods; 2011–2045, 2046–2079, and 2080–2112. In regard to the findings, the upward trend of the potential evapotranspiration in parallel with the downward trend of the precipitation in the next 102 years in three scenarios to the base timescale was transparent. The frequency of the SPEI in the base month indicated that 17.02% of the studied months faced the drought. Considering the scenarios of climate change for three 34-year periods (i.e., 2011–2045, 2046–2079, and 2080–2112) the average percentages of potential drought occurrences for all the stations in the next three periods will be 8.89, 16.58, and 27.27 respectively under the B1 scenario. While the predicted values under the A1B scenario are 7.63, 12.66, and 35.08%respectively. The relevant findings under the A2 scenario are 6.73, 10.16, 40.8%. As a consequence, water shortage would be more serious in the third period of study under all three scenarios. The percentage of drought occurrence in the future years under the A2, B1, and A1B will be 19.23%, 17.74%, and 18.84%, respectively which confirms the worst condition under the A2 scenario. For all stations, the number of months with moderate drought was substantially more than severe and extreme droughts. Considering the A2 scenario as a high emission scenario, the analysis of SPEI frequency illustrated that the proportion of dry periods in regions with humid and cool climate is more than hot and warm climates; however, the duration of dry periods in warmer climates is longer than colder climates. Moreover, the temporal distribution of precipitation and potential evapotranspiration indicated that in a large number of stations, there is a significant difference between them in the middle months of the year, which justifies the importance of prudent water management in warm months.

duration of dry periods in warmer climates is longer than colder climates.Moreover, the temporal distribution of precipitation and potential evapotranspiration indicated that in a large number of stations, there is a significant difference between them in the middle months of the year, which justifies the importance of prudent water management in warm months.

Introduction
Drought as one of the natural disasters has enormous socio-economic consequences all over the world.Its long-lasting effects on agriculture, the ecosystem, and the economy over broad regions have been reported in several investigations [1,2].Since determining the exact dates of the beginnings and endings of droughts is never possible, estimating the effects and losses of drought is difficult.Furthermore, droughts could cause or accentuate other dangerous events such as floods, storms, and hurricanes [3].This expensive phenomenon is driven by a variety of unpredictable factors which in turn enhances the complexity of this poorly understood phenomenon [4,5].On the other hand, climate change as a reaction of nature to human activities is projected will likely increase in the future due to the global warming effect, which deteriorates water shortage and threatens human lives and livelihoods [6].Accordingly, climate change influences the ecosystem in terms of changing precipitation patterns and air temperatures, which may cause acute droughts.In the presented work [7], climate change would intensify natural disasters such as drought or flood.These consequences are visible in either shortterm or long-term processes [8,9].
Given that, identifying, and discovering changes in drought and finding their culprits are complicated.Moreover, high spatial and temporal variability of drought conditions results in uncertainty in forecasting drought.These circumstances necessitate the study of this natural phenomenon in every area which has been done in several investigations [10].The climate change effect evaluation conducted on meteorological drought indices including Standard precipitation index (SPI) and Standardized precipitation-evapotranspiration index (SPEI) in Varamin plain [11].Kamali and Asghari (2023) investigated the impact of future drought events under climate change conditions on groundwater storage in Najafabad aquifer in Iran [12].Muthuvel et al. (2023) assessed global concurrent drought traits and their effects on maize yield under climate change in the USA, China, Brazil, Mexico, Argentina, France, and India [13].The impacts of climate warming on the future hydrological droughts employed in 6,688 catchments by Gu et al. [14], and revealed that precipitation stress controls the development of historical droughts over *87% of global catchments.Predictability of meteorological drought in the semi-arid region based on the standardized precipitation index (SPI) was investigated by Elbeltagi, et al. [15].An evaluation established for the future drought in the west of Iran by SPI under the A1B, B1 and A2 emission scenario of HadCM3 for the period 2011-2030 [10].The result of their study indicated that the severity of drought will increase with longer return period.Quantifying the effect of climate change on characteristics of drought in Australia, Brazil, China, Ethiopia, India, Spain, Portugal and the USA showed that adaptation to drought risk will be vital in these regions for the first half of the twenty-first century [16].Assessing features of future drought conditions over South Korea by [17] demonstrated an increase of 6% in frequency and severity of drought during the dry season.In the presented work [18], the authors investigated drought risk under climate and land use changes in the Frome catchment in the UK.In this survey, the drought indices predicted an increase in the severity of future drought events under the high emission scenarios.The study of [19] revealed an increasing and worsening drought event in the future than the last three decades for Northern Ethiopia.Monitoring drought patterns in the future climate highlighted the urgency for breeding drought-tolerant lines to secure crop productivity in rainfed regions such as Australia.Application of potential evapotranspiration (PET) in assessing influences of climate change on drought variation in the Andean high mountain basin illustrated that the evaluation of droughts in complex terrain in tropical regions should be promoted, where the representation of convection is the main limitation of global climate models [20].In other work [15], the authors projected SPI and SPEI under adjusted outputs of global climate models in Turkey.The projections indicated fewer drought events for the near-future period of 2016-2040, with no potential extreme drought events.[21] compared indices of SPEI, SPI and Streamflow Drought Index (SDI) in identifying future drought characteristics in Cheongmicheon watershed, Korea, under general circulation models.According to the obtained results, the SPEI showed the most robust pattern with significant increases in both drought frequency and severity.
In regard to the previous investigations [22,23], it can be inferred that drought characteristics alter over time due to climate changes and among regions with different climate conditions [24].Accordingly, forecasting and monitoring drought would be worthwhile and beneficial in management of this phenomenon and would reduce losses due to it.Although several studies have addressed similar topics, it should be mentioned that the vast majority of them are based on regional scales for a short period of time.Whereas, the investigation of drought for a long time in large scales considering different climatic conditions under distinct climate change scenarios has been paid less attention.This kind of study is more benefit for an area with acute water shortage such as Iran which faces long-lasting droughts for several times.Given that, the present investigation could have an advantage in terms of considering distinct climates over a large scale.Moreover, evaluating drought using the SPEI for a long term in different climates can justify the importance of the present study.The findings of this investigation could be used in managing water resources and administering sufficient policies for tackling the future water shortage crisis.

Material and methods
Geographically, Iran is located at 25.05 to 39.78 Northern latitude and 44.08 to 63.30 Eastern longitudes.This study relied on the records of daily dataset collected from 30 synoptic stations for a 20 years period (1992-2010) for the assessment of the impacts of changes in climate conditions on drought.These stations were mainly considered for the quality of the data, as well as the uniformity in the period of data collection.Although some stations have data for a longer period of time, considering stations with the same time scale and no data missing matters to comparing the results under a reliable condition.There is no hesitation that any filling data technique has an error and will affect the findings, accordingly, to avoid any error related to missing data, this period of time was considered which had no data gap.In other words, datasets based on poor gauge coverage can produce substantial uncertainty when gaps are filled with data from different sources based on some climatology statistics [25].Another consideration is that the stations are in different climatic locations which give a good sense of regional spatial distribution.Fig 1 depicts the distribution and location of the considered stations, as well as their climatic predispositions as per the Koppen system of climate classification [26].
Northern Iran is mainly the farming area of the country to the favourable weather condition in the area.This necessitated the citation of several weather stations in northern Iran than the other areas of the country.The assessment of drought can be carried out by a variety of indices, however accuracy of each index changes in different regions.For instant, high accuracy of the SPEI compared to other indices in different provinces of Iran has been reported by Zarei et al. [27].The benefit of SPEI has been indicated in several investigations as well.Uddin et al (2020) mentioned that SPEI performs better than SPI because evaporative demand has a positive impact on defining drought conditions [28].Wang et al. (2021) revealed that SPEI can be used to delineate spatial-temporal evolution of drought, drought characteristics, and impacts of drought at the regional and global scales which confirms the advantages of this drought index [29].Accordingly, the SPEI was used in this study to evaluate drought conditions.
The determination of drought threats using the SPEI method was first introduced by [7] as a method that depends on two major parameters-temperature and precipitation.Being that temperature is the major variable in PET calculation, drought conditions are mainly determined by considering the variation between precipitation (P) and PET: Estimating accumulated probabilities of D i by probability density function is essential to calculate SPEI [30].Due to negative values in D i , the probability density function of a three- parameter log-logistic is needed to calculate the SPEI as the following [7]: where α, β and γ are scale, shape, and origin parameters, respectively, for D values (γ�D�1).
The probability distribution function of the D series, according to the log-logistic distribution, is given by: The F(x), values for the D series at different time scales adapt very well to the empirical F(x) values at the different observatories, independent to the climate characteristics and the time scale of the analysis.
With F(x), the SPEI can easily be obtained as the standardized values of F(x).For example, following the classical approximation of [31]: where; W ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi À 2 lnðpÞ p for p � 0:5 P presents the computed D value exceeding probability.P = 1-F(x).if >0.5, then P is substituted by (1-p) and the sign of the resultant SPEI is reversed.
The SPEI mean and standard deviation values are zero and one, respectively.SPEI is an indication for standardized variable and thus it is comparable with other SPEI magnitudes over time and space.An SPEI of 0 signifies a corresponding value to 50% of the cumulative probability of D, this is according to a log-logistic distribution [32].In the present investigation, the SPEI was run on package of the R 3.63 software.
Since the definition and conception of SPI is applied in the development of SPEI, different classes of drought intensity are based on recommended values by [33].Table 1 reveals the classification of the drought categories with respect to SPEI.
For describing and analyzing the drought characteristics, the Run theory was used.In this theory, severity, duration, and intensity are considered as three important characteristics of drought.Fig 2 depicts these traits with X 0 as threshold of X t [34].
First, the climatic variables were computed by LARS-WG for the future period as this is needed for the calculation of future PET.LARS-WG was calibrated using maximum and minimum temperatures (Tmax and Tmin) and precipitation during base period of 1992-2012 for each station.By this way, these three variables (Tmax, Tmin and precipitation) were simulated for the future periods.For the calibration of LARS-WG, two comparison statistics, mean bias error (MBE) and normal root mean square error (NRMSE) were applied, which are defined as [35]: where, X are the observed values and Y are the simulated values, � Xand � Y indicate the means of X and Y values respectively, n is the number of data.
After a good calibration of the developed model for the base period, the mentioned climatic variables were predicted for the future periods considering three scenarios, A1B, A2 and B1.
The assumption made by the B1 scenario relies on an endurable world, that is, quick changes in economic developments, advancement of human rights equality and a consideration to environmental protection.This assumption considers that greenhouse gas emissions can be controlled and a pollutant controlling programme will be in place for the industries.The assumption made by the scenario of A1B relies on a wealthy world with a quick economic development (about 3percent per year), decrement in population (about 27% per year), quick technological development, convergence of culture and economy and a principle decrease in regional differences.The assumption made by the scenario of A2 relies on the assumption of a separate world.Different cultural characteristics in various areas of the world are expected to cause differences in the change of climate parameters [36].Hence, the climatic variables were forecasted for three different periods (each having 34 years), 2011-2045, 2046-2079, and 2080-2112 under various scenarios.
In the second step, the amount of PET was calculated by the Hargreaves-Samani (HS) method.Due to the fact that the SPEI considers both precipitation and temperature, thus applying a method based on temperature can reflect the impact of temperature, which is completely transparent in the HS.Moreover, through comparing influence of climatic parameters on describing drought in Iran reported that temperature and precipitation parameters had the most influential effect on the fluctuation of drought severity in Iran [37].Furthermore, previous works suggested this method in Iran [36,38].Therefore, PET was estimated by the HS equation which utilizes maximum and minimum temperatures, mean temperature and solar radiation as follows: where, PET is the potential evapotranspiration (mm), T is the mean air temperature (˚C), TR is the difference between Tmax and Tmin (˚C), and Ra is the extraterrestrial radiation (MJm - 2 day -1 ).The Ref-ET software was implemented for estimating PET.Finally, considering the precipitation parameter during future years, the magnitude of the SPEIs were predicted for the future 102 years from 2011 to 2112 under each scenario.
For assessing the changes in the trend of SPEI between 2011 and 2112, the Mann-Kendall test was used.The Mann-Kendall test is applied to see whether Y values tend to decrease or increase with T (monotonic change).It can be hypothesized as the concentration independent of time of the chemical parameters [39] or H 0 : no correlation between x and y; Ha: x and y are correlated.To apply the test, Kendall's S statistic is calculated from the Y, T data pairs: where, P = the number of Y i <Y j for all i<j, M = the number of Y i >Y j for i<j: ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi

Modeling results and analysis
Evaluating the accuracy of LARS-WG performance in simulating meteorological variables plays a key role in predicting weather parameters for the future periods.For this purpose, the consequences of statistical error indices in the base period are presented in Table 2.
According to the information of Table 2, it can be inferred that the calibration of the model was successful with high accuracy.The values of NRMSE and MBE confirm the lowest error in simulating weather variables.Although there are significant differences among the values of NRMSE and MBE in few stations, in general, the results of Table 2 verify proper performance of the model in creating connection between observed and simulated values.
Then, the prediction of the weather parameters for the upcoming period was done as earlier described in the materials and methods section.PET values were calculated using Eq 8 while the monthly spatial precipitation and PET distribution for the considered period and for up to Based on the assumptions of the IDW interpolation, features are indistinguishable if they are placed close to each other compared to the ones farther apart.The IDW strategy has been suggested to be used for data with a short range of variation [36].It is obvious from Comparing the future PET under the scenarios indicate that the trend of spatial changes of PET under the A2 and A1B scenarios are similar to the PET in the base period, whereas there is a narrow margin between the PET changes under the B1 scenario and the base period.Nevertheless, the increasing trend of the PET in the three scenarios to the base timescale is transparent.These findings are reported by [36] that can confirm the current consequences.However, the spatial distribution of the precipitation pattern is completely the reverse of the PET, which is reasonably due to the identity of both the parameters.Contrary to the PET values, the lowest precipitation falls in the south-eastern parts of the country.This is noticeable that the category related to the highest precipitation under the A2 scenario is significantly less than other scenarios.In other words, the vast area of the country will receive precipitation less than 23.6 mm per month.However, this condition under the A1B and B1 scenarios seem better.Reduction of future precipitation in Iran has been reported by [40][41][42], that are in line with the current results.Table 3 depicts the drought characteristics of each station in the base period.Describing drought was based on the recommendation of Edwards and McKee (1997) given in Table 1.The information of this table shows the percentage of dry and wet periods in the base period.For instance, value 1.32 percent in Ahvaz refers to the fact that the SPEI in 3 months out of 228 months (19 years) are less than -2, which refers to an extreme drought (Table 1).Overall, 17.02 percent of the studied months in the base period faced the drought, whereas 17.84 percent were in the wet condition.The most intensive drought was observed in Boshehr, by 439 percent.The highest number of dry months belonged to Birjand with 25.44 percent.It is noticeable that regions with high precipitation like Rasht or Gorgan have also experienced drought periods.The similar results are reported by [37], that support the present findings.
For presenting the calculated SPEI values for the future period, the percentage of months that SPEI values are within each category to all the considered months (e.g., 102 years or 1224 months) was counted.The findings under each scenario are shown in Tables 4, 5, and 6.In these tables, the number of months in each category has been presented as a percentage for each studied station.
From Table 4, it can be seen that the response of each station to future climate change is completely different during the studied period.Nevertheless, for most of the stations, the percentage of drought events in the third period (2079-2112) is more than the other periods.The average percentages of potential drought occurrences for all the stations in the next three periods are 8.89, 16.58, and 27.27 respectively.These values refer that the probability of drought will sharply increase over time.On average, about 17.74 percent of months over all the stations are likely to experience drought throughout the future years up to 2112.This finding reveals that the frequency of drought events will increase under the B1 scenario compared to the base period.Among them, the highest percentage of dry months is observed in Shiraz with 19.69%, while the lowest proportion is related to Bandarabas with 14.79 percent.Comparison of dry months demonstrates that the great number of months with extreme drought is seen in Yazd with 1.23 percent.For all stations, the number of months with moderate drought is substantially more than severe and extreme droughts.In contrast, the consequences of the humid events indicate that approximately 18.55 percent of months in all stations will experience wet months.Esfahan will experience the highest frequency of humid months over the upcoming period; however, the number of months with extreme wet is seen in Boshehr with 2.37%(29 months).It is remarkable that for all stations, the occurrence probability of moderate humid months is more than that of severe and extreme humid months.Although, the number of months with extreme wet or extreme drought is substantially less than others, it should be noted that losses and damages of extreme wet or drought would be considerably more than    other situations.Thus, crisis management in each region should be designed for extreme conditions.
From the data given in Table 5 it can be inferred that like the results of Table 4, for all the regions the frequency of drought in the third period is substantially more than the previous periods.This value is about 35.08 percent, which is higher than the B1 scenario, by a broad margin.As an overall trend, the results of Table 5 indicate that, the frequency of dry months in Iran over the studied period (from 2011 to 2112) will be about 18.46 percent, which is more than the observations in the B1 scenario (e.g., 17.74 percent).Thus, it can be concluded that the proportion of predicted drought under the A1B scenario is more than the B1 scenario.Subsequently, this value is far more than the number of droughts in the base period.This finding is in agreement with the survey of [43], who reported that the emission scenario of B1 always produces less severe extreme disasters.This is because the emission scenario of B1 assumes that the future world will be sustainable with quick variations in economic structures and more efforts in introducing clean technologies.
Detailed consequences illustrate that the number of months that will be affected by moderate drought is considerably more than others.Arak can be introduced as a station that will face the highest number of moderate droughts over the next period, with 199 cases (16.26%).However, the most proportion of months with extreme drought will be seen in Ardebil with 1.31percent, which can be a serious problem for water resources of this region.Overall, the highest frequency of dry months will belong to Shiraz, 259 months from 1224 months (21.16%).Thus, it can be mentioned that the probability of drought incidence in Shiraz would be the most among all the stations.
Evaluating the results related to humid periods depict that for all stations, the number of months with moderate wet is significantly more than those with severe and extreme wet.Moreover, the number of months with extreme wet is considerably less than the severe wet.Looking in more detail, it can be seen that the highest frequency of humid months will be observed in Ghazvin with 20.10 percent, whereas the lowest frequency would be related to Hamedan, with 16.01percent.The findings of wet intensity show that the greatest number of extreme wets will occur in Boshehr, with 2.46 percent in comparison to 0.90% for Ghazvin as the lowest number.
The information given in Table 6 depicts that the frequency of droughts in the third period are dramatically more than the previous periods.This difference is completely significant.This amount is about 40.81 percent in comparison to 6.72 and 10.16percent for the first and second periods respectively.
It can be seen from Table 6 that the number of dry months under the A2 scenario (19.23percent) is more than the B1 and A1B scenarios (17.74 and 18.46percent).Due to the fact that the A2 scenario is introduced as a high emission scenario which describes a very heterogeneous world, this result seems to be reasonable and rational.Furthermore, this finding is consistent with other studies.For instance, in the survey of Akbari et al. (2016) [44], the highest temperature, and the lowest precipitation prognosticated belonged to the A2 scenario.Furthermore, this result corresponds to the results of [45].They indicated that a drought event with more intensity, duration and frequency than the base period can occur in future periods under the A2 scenario in northwest of Iran.According to the study of Modarres et al. (2018) [46], the average of the rainfall events will reduce in evaluating future extreme rainfall changes in the north of Iran as well.It can be an alarm in applying effective strategies for the future water resources in this area.An increasing trend of drought under the A2 scenario has been declared by [47], which is in line with the present forecasted consequences.The conducted research by [48] for assessing the future temperature in Iran, researchers declared that the maximum temperature changes will be increased more severely based on the A2 scenario compared to other scenarios, which is in line with the current investigation.Although in some investigations like [21], an increased precipitation under the B1 scenario compared to other scenarios has been reported, it should be mentioned that the time scale is different in this study.Moreover, since the changes of both precipitation and PET are taken into account for calculating SPEI, the interpretation of SPEI can be distinct from precipitation.
In parallel with results of the B1 and A1B scenarios, the most frequency of dry months is observed in Shiraz with 22.63percent, while this value for the B1 and A1B scenarios are 19.70 and 21.16 percent (Tables 4 and 5).1.23percent of months in Bandarabbas expected to experience extreme droughts, which can be labeled as the driest station throughout the next years.This result matches with the survey of [37].Due to the proportion of dry months, it can be inferred that approximately over 16.42 percent of months of each station will be experiencing droughts during the next years, which is alarming for Iran.Although the number of months with extreme droughts is much more than other categories, it should be declared that losses and damages of extreme intensity droughts are dramatically more destructive than others.In other words, a short-time drought period with extreme intensity can have long-term repercussions.
The highest figure of humid months under the A2 scenario belongs to Ghazvin with 19.93percent, which is similar to the results of the A1B scenario.It is prominent that for all the stations, the number of months with moderate wet is significantly more than the cumulative amount of extreme and severe wets.The highest number of months with extreme wet is related to Bosheher with 2.30 percent.
For examining the changes in the trend of SPEI over the studied period from 2011 to 2112, the results of Mann-Kendall test (Z) are provided in Fig 4 .The aim of this test is to statistically evaluate if there is a monotonic increase or decrease trend of the studied variables over the time.It is noticeable that Z�|1.96| refers to the significance level 0.05, while Z�|2.58|implies to the significance level 0.01.significantly.So that, for most of the stations, this trend is significant in the level of 0.01.The highest downward trend for a large number of stations is observed in the A2 scenario.Thus, it can be declared that the amount of SPEI will drop sharply over time and this decrease under the A2 scenario will be the most intensive.
As a consequence, it can be inferred that the A2 scenario predicts the most serious drought in comparison to other scenarios.Moreover, since the survey of [36] demonstrated that the certainty of the predicted PET values under the A2 scenario is higher than the other scenarios, the results of the A2 scenario were examined for extracting the characteristics of SPEI based on the Run theory,.Findings are provided in Table 7 for each station for the next 102 years.
As an overall trend, comparing the number of dry periods for all the stations indicate that the proportion of dry periods in humid and cold regions is more than warm climate.For instance, 288 dry months are expected for Gorgan, in comparison to 218 dry months for Zahedan.However, in regard to long dry periods, the longest dry periods are expected to be observed in warmer stations; 8 in Gorgan compared to 25 in Zahedan.The maximum deficit can also be seen in warm climates such as Zahedan, with 27.36.This consequence can be seen It should be mentioned that uncertainty in the prediction of future events is inevitable, it is significantly affected by the type of data used in terms of quality and quantity aspects, and also the structure of downscaling model [36,49].Regarding the data, using dataset with no missing data will reduce uncertainty, which was done in this study.Apart from that, the SPEI index assimilates more factors to quantify drought characteristics which affects less uncertainty [25].From aspect of model, considering the calibration results, it can be claimed that the present model could discover connections of weather parameters with high accuracy, thus the performance of the model is acceptable.Moreover, different scenarios provided a chance to consider different assumptions in prediction which was compared.There is no hesitation that removal of uncertainty in modelling is impossible, nevertheless, the above-mentioned reasons can confirm reliability of the findings in this study.
For more understanding, the monthly changes of PET and precipitation predicted under the A2 scenario for each station are provided in Fig 5.It seems to be crystal clear that for the vast majority of stations, the area of PET is significantly more than the area of precipitation.The most difference can be seen in the middle months of the year (e.g., in warm months), which is a logical and rational consequence.These findings are consistent with the report of [50], who concluded that the magnitudes of PET in Iran increases from January to July, then drop from August to December, and all stations reach their maximum levels in July.Since the highest temperature takes place in warm months, subsequently the amount of PET would increase as warming accelerates.However, the probability of precipitation occurrence would drop.These findings can be observed in regions with hot climates such as Semnan, Yazd, Zahedan, and Kerman which are located in the southern part of Iran.[37] also declared that the southern regions of Iran are more severely affected by drought.These results refer to the importance of suitable water management in warm months and especially preserving and harvesting precipitation during the months with high precipitation.Contrary to these regions, the amount of precipitation (110.80 mm/month) in Rasht is approximately 26 percent more than the PET (87.51 mm/month).These consequences are in line with the findings of [51].
In general, it can be concluded that water shortage in warm climates would be much more serious than in other climates.Despite insufficient precipitation being the principal cause of water scarcity and drought, non-meteorological drivers such as rapid population growth, poor water resource management, and unsuitable cropping patterns also play critical roles [51].

Conclusion
Evaluating the amount of precipitation and PET for the base period and the next 102 years depicted that the highest PET and the lowest precipitation are related to the southeast of the studied area.Prediction of drought by SPEI under climate change scenarios indicated that the vast majority of stations will be involved in drought events.However, regions with hotter climates will experience more intensive drought.Comparing forecasted drought under the A2, A1B and B1 scenarios showed that 19.23% of months over the next 102 years under the A2 scenario will struggle with water shortage and dry periods, which is more than other scenarios and the base time scale.The study highlights these findings: the number of dry periods in cold and cool terrains would be more than hot climates, whereas their intensity would be less than hot and warm climates.Furthermore, the temporal distribution of PET and precipitation revealed that for most of the stations, the precipitation was considerably less than the PET, especially in the middle months of the year.This situation in the southern regions is likely to be more significant.As a consequence, the severity and number of future droughts are expected to deteriorate particularly in warmer climates in terms of both meteorological and hydrological aspects, which highlights the importance of current water resources management for transmitting these natural resources to the next generations.According to the findings, this study broadens the knowledge on the future drought condition over Iran with acute water shortage crisis.Based on the observation, different regions will encounter distinct drought circumstances in the future which necessitates implementing specific plans to address water deficit in each area.

*
Positive values of Z show upward trends while negative values of Z indicate downward trends.
the next 102 years (2011 to 2112) under different climate scenarios is shown in Fig 3.These figures were produced with ArcGIS 10.3 using the inverse distance weighted (IDW) procedure.
Fig 3 that the Southern part of Iran has higher PET values (ranges from 124.1 to 145.6 mm/month) than the northern part (ranges from 82.47 to 104 mm/month).

Fig 3 .
Fig 3. Spatial distribution of the monthly PET and precipitation in the base period and the future under the A1B, A2 and B1 scenarios.https://doi.org/10.1371/journal.pone.0290698.g003

Fig 4
indicates that there is a downward trend in the SPEI changes in all the stations, except Esfahan under the B1 scenario.In other words, the value of SPEI over time will decrease

Fig 4 .
Fig 4. The result of Mann-Kendall test on SPEI trend over the whole period.https://doi.org/10.1371/journal.pone.0290698.g004 Fig 5 presents the area charts of PET and precipitation in a monthly time scale for the next 102 years.

Table 2 . The results of LARS-WG test in prediction of weather variables after calibration in the base period.
https://doi.org/10.1371/journal.pone.0290698.t002

Table 7 . Characteristics of dry and wet periods based on SPEI under the A2 scenario. Station Dry period Wet period Number of periods Longest period Maximum deficit Total deficit Number of periods Longest period Maximum surplus Total surplus
other stations such as Ghom with a hot climate and Tabriz with a cold climate.In regard to wet periods, the results are the same, i.e. the highest frequency of wet periods is related to cold climates, whereas the lowest belongs to warm climates.However, the duration of dry periods in warm climates is longer than that of cold climates.The information given for Rasht and Gorgan in comparison with Zahedan and Ahvaz can confirm these declarations.Sobhan et al. (2019) also reported that the southern stations of Iran are more exposed to drought than others. in